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The theoretical predictions of the primordial abundances of elements in the big-bang nucleosyn- 
thesis (BBN) are dominated by uncertainties in the input nuclear reaction rates. We investigate the 
effect of modifying these reaction rates on light element abundance yields in BBN by replacing the 
thirty-five reaction rates out of the existing eighty-eight. We have studied these yields as functions 
of evolution time or temperature. We find that using these new reaction rates results in only a little 
increase in helium mass fraction over that obtained previously in BBN calculations. This allows 
insights into the role of the nuclear reaction rates in the setting of the neutron-to-proton ratio during 
the BBN epoch. We observe that even with considerable nuclear physics uncertainties, most of these 
nuclear reactions have minimal effect on the standard BBN abundance yields of 6 Li and 7 Li. 
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I. INTRODUCTION 

Primordial nucleosynthesis took place just a few mo- 
ments after the big-bang [l| and the universe evolved 
very rapidly allowing only the synthesis of the lightest 
nuclides such as D, 3,4 He and 6 - 7 Li. In addition to these 
stable nuclei some unstable, or radioactive, isotopes like 
tritium or 3 H and 7 ' 8 Be were also produced during the 
big-bang nucleosynthesis (BBN) . These unstable isotopes 
either decayed or fused with other nuclei to make one of 
the stable isotopes. It lasted for only about seventeen 
minutes (during the period from three to about twenty 
minutes from the beginning of space expansion) and af- 
ter that, the temperature and density of the universe 
fell below that which is required for nuclear fusion and 
prevented elements heavier than beryllium from form- 
ing while at the same time allowing unburned light ele- 
ments, such as deuterium, to exist. The abundances of 
these nuclides are probes of the conditions of the universe 
during the very early stages of its evolution. The condi- 
tions during the BBN are believed to be well described 
in terms of standard models of cosmology, nuclear and 
particle physics, which determine the values of tempera- 
ture, nucleon density, expansion rate, neutrino content, 
neutrino-antineutrino asymmetry, etc. Deviations from 
the BBN test the parameters of these models and con- 
strain nonstandard physics or cosmology that may alter 
the conditions during BBN Sensitivity to several 

parameters and physics input in the BBN model have 
been investigated thoroughly in the past [5MT0I] . 

The nuclear reaction rates < av > in reaction network 
calculations, where a is the nuclear-fusion cross section 
and v is the relative velocity between the participant nu- 
clides, are the most important inputs for modeling the 
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BBN and stellar evolution. Although v is well described 
by a Maxwell- Boltzmann velocity distribution for a given 
temperature T, the cross section a can only be obtained 
from laboratory experiments, some of which are not as 
well known as desired • Several factors influence the 
measured values of the cross sections and the theoretical 
estimates of the thermonuclear reaction rates depend on 
the various approximations used. In the network calcula- 
tions one needs to account for the Maxwellian-averaged 
thermonuclear reaction rates and difference [lj]> Ell m 
these rates affects the description of elemental synthe- 
sis in the BBN or in stellar evolution. In this work, we 
consider the effect the nuclear reaction rates on the pri- 
mordial abundances of elements. 



II. BIG-BANG NUCLEOSYNTHESIS 
REACTION NETWORK 

The predictions of the standard BBN theory depend 
on low energy nuclear cross sections and on three addi- 
tional parameters, the number of flavours of light neu- 
trinos (N u ), the neutron lifetime (r„) and the ratio of 
baryons to photons (77 = ns/n 7 ) in the universe (l3l.[l4j]. 
In its standard N v = 3.0 form, primordial nucleosyn- 
thesis is, therefore, a one parameter theory, depending 
only on r\ as t„ = 885.7(8)s is an accurately measured 
quantity. One of the most essential predictions of the 
standard big-bang model is the synthesis of light ele- 
ments in the primordial universe. There are three pillars 
of big-bang cosmology which are the Hubble expansion, 
cosmic background radiation (CBR) and the BBN. The 
BBN examines back to earlier times of the universe than 
the other two and it also has to deal with nuclear and 
particle physics along with cosmology. Although Hub- 
ble expansion can also be used by some other alternative 
cosmological theories (e.g. steady state), the evidences of 
CBR and BBN observations lead cosmologists to a uni- 
verse which was very hot and dense in the very beginning. 



2 



The Friedmann-Robertson- Walker cosmological model is 
the standard scenario underlying the BBN theory. More- 
over, the solution of the Einstein equations lead to an 
isotropic and homogeneous universe, so the uniformity 
of the CBR temperature, which is T = 2.7277 ± 0.002°K 
across the sky, as well as the success of the standard BBN 
theory serve to validate this approximation. The primor- 
dial harvest of light elements is determined by the time 
during the expansion of the universe. It is possible to 
characterize the BBN in general, that the paradigm most 
frequently uses the Friedmann equation to relate the big- 
bang expansion rate, H, to the thermal properties of the 
particles present at that epoch. During expansion, the 
rates of the weak interactions that transform neutrons 
and protons, and the rates of the nuclear reactions that 
build up the complex nuclei are involved. 

The reactions which happened at the duration of BBN 
can be organised into two groups, firstly, the reactions 
that interconvert neutrons and protons which are n + 
e + -R> p + P e , p + e~ n + v e and n p + e~ + v e and 
secondly, the rest of the reactions. The first group can 
be expressed in terms of the mean neutron lifetime and 
the second group is determined by many different nuclear 
cross section measurements. The formation of deuterium 
begins in the process ofp + n-s->£> + 7. This reaction is 
exothermic with an energy difference of 2.2246 MeV, but 
since photons are 10 9 times more numerous than protons, 
the reaction does not proceed until the temperature of 
the expanding universe falls to about 0.3 MeV, when the 
photo-destruction rate is lower than the production rate 
of deuterons. When the deuteron formation starts some 
further reactions proceed to make 4 He nuclei. 



3 H + p -> 4 He + 7 
D +p -> 3 He + 7 
3 He + 4 He + 7 



(1) 



Both light helium 3 He and normal helium 4 He are 
formed along with the 3 H . Since helium nucleus binding 
energy is 28.3 MeV and more bound than the deuterons 
and the temperature has already fallen to 0.1 MeV, these 
reactions can be photoreactions and only go one way. 
The following four reactions also produce 3 He and 4 He 
and they usually go faster since they do not involve the 
relatively slow process of photon emission. 



D + D- 
D + D 
D - 

D - 



3 He 
3 H 
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Eventually the temperature gets so low that the electro- 
static repulsion of the deuterons and other charged parti- 
cles causes the reactions to stop. The deuteromproton ra- 
tio when the reactions stop is quite small, and essentially 



inversely proportional to the total density of protons and 
neutrons (to be precise, goes like the -1.6 power of the 
density). Almost all the neutrons in the universe end up 
in normal helium nuclei. For a neutron:proton ratio of 1:7 
at the time of deuteron formation, 25% of the mass ends 
up in helium. Deuterium peaks around 100 seconds after 
the big-bang, and is then rapidly swept up into helium 
nuclei. A very few helium nuclei combine into heavier 
nuclei giving a small abundance of 7 Li coming from the 
big-bang. 3 H decays into 3 He with a 12 year half-life so 
no 3 H survives to the present, and 7 Be decays into 7 Li 
with a 53 day half-life and also does not survive. The 
uncertainties for the reactions 3 He + 4 He — > 7 Be + 7, 
3 H + 4 He -^ 7 Li + -f and p + 7 Li -> 4 He + 4 He may lead 
to about 50% uncertainty in the predicted yield of 7 Li. 
In the present work, we replace thirty-five Maxwellian- 
averaged thermonuclear reaction rates by new ones in the 
Kawano/ Wagoner BBN code [l5[ and study its effect on 
the primordial abundances of elements. 



III. THERMONUCLEAR REACTION RATES 

The twelve most important nuclear reactions which 
affect the predictions of the abundances of the light 
elements [ 4 He, D, 3 He, 7 Li] are n— decay, p(n,"f)d, 
d{p,-/) 3 He, d(d,n) 3 He, d(d,p)t, 3 He(n,p)t, t(d,n) 4 He, 
3 He(d,p) 4 He, 3 He{a,j) 7 Be, t(a,j) 7 Li, 7 Be(n,p) 7 Li 
and Li{p, a) He. Instead of cross sections a, the nuclear 
reaction inputs to BBN take the form of thermal rates. 
These rates are computed by averaging nuclear reaction 
cross sections over a Maxwell-Boltzmann distribution of 
energies. The Maxwellian-avcragcd thermonuclear reac- 
tion rate < av > at some temperature T, is given by the 
following integral 161 ]: 



(3) 

where E is the centre-of-mass energy, v is the relative ve- 
locity and [i is the reduced mass of the reactants. At low 
energies (far below Coulomb barrier) where the classical 
turning point is much larger than the nuclear radius, bar- 
rier penetrability can be approximated by exp(— 2-7r£) so 
that the charge induced cross section can be decomposed 
into 



a(E) 



S(E) exp(-27<) 
E 



(4) 



where S(E) is the astrophysical S'-factor and £ is the 
Sommcrfcld parameter, defined by 



c 



ZiZ 2 e 2 

h v 



(5) 



where Z\ and Z 2 are the charges of the reacting nuclei 
in units of elementary charge e. Except for narrow reso- 
nances, the 5-factor S(E) is a smooth function of energy, 
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which is convenient for extrapolating measured cross sec- 
tions down to astrophysical energies. In the case of a nar- 
row resonance, the resonant cross section a r {E) is gener- 
ally approximated by a Breit-Wigner expression. 

TABLE I: The nuclear reactions with modified thermonuclear 
reaction rates for BBN are tabulated along with the applicable 
Tg (in units of 10 9 °K) ranges and references. Rows without 
Tg ranges show applicability for the entire range. 





Tg < 0.8 


[lJJ 


r 9 


> 0.8 [12] 


d(d,n) 3 He 


Tg < 3.0 


[12] 


r 9 


> 3.0 [12] 


d{d,p)t 


T 9 < 3.0 


[UJ 


Tg 


> 3.0 [12] 


d(a,j) 6 Li 


[12] 








t{d,n) 4 He 


Tg < 0.5 


[17] 


T, 


> 0.5 [12] 


t(a, 7) 7 Li 


Tg < 8.0 


[UJ 


Tg 


> 8.0 [12] 


3 He(n,p)t 


Tg < 3.0 


[17] 


Tg 


> 3.0 [18] 


3 He(d,p) 4 He 


Tg < 2.0 


[17] 


T 9 


> 2.0 [18] 


3 He( 3 He,2p) 4 He 


[12] 








3 He{a, 1 ) 7 Be 


Tg < 8.0 


[UJ 


Tg 


> 8.0 [12] 


4 He(an,-f) 9 Be 


[12] 








A He{aa,-y) 12 C 


[12] 








6 Li{p,-/) 7 Be 


[12] 








6 Li(p,a) 3 He 


[12] 








7 Li(p,a) 4 He 


T 9 < 7.0 


[UJ 


T 9 


> 7.0 [12] 


7 Li(a,~ t ) 11 B 


[12] 








7 Be(n,p) 7 Li 


T 9 < 0.2 


[UJ 


Tg 


> 0.2 [18] 


7 Be{ Pn fB 


[12] 








7 Be(a, 7 ) 11 C 


[12] 








9 Be(p,~f) 10 B 


[12] 








9 Be(p,da) i He 


[12] 








9 Be(p,a) 6 Li 


[12] 








9 Be(a,n) 12 C 


[12] 








10 B(p,7) n C 


[12] 








10 B{p,a) 7 Be 


[12] 








n B(p,7) 12 C 


[12] 








lx B(p,aa) A He 


[12] 








12 C(p, 7 ) 13 iV 


[12] 








12 C(a, 7 ) 16 


[12] 








13 C(p,y) 14 N 


[12] 








13 C(a,n) 16 


[12] 








13 iV(p,7) 14 


[12] 








14 iV(p, 7 ) 15 


[12] 








15 N(p r/ ) la O 


[12] 








15 N(p,a) 12 C 


[12] 









The neutron induced reaction cross sections at low en- 
ergies can be written as [19( 

<e) = tm. ( 6 ) 

V 

facilitating extrapolation of the measured cross sections 
down to astrophysical energies, where R(E) is a slowly 



varying function of energy [201 ] and is similar to an S- 
factor. 

The compilations concerning specifically standard 
BBN reaction rates were performed by Caughlan et al. 
[llj and Smith et al. [Hj]. We are using rates for some 
BBN reactions from the newest compilations by Angulo 
et al. [IH and Descouvemont et al. [TtJ which are meant 
to supersede the earlier compilations. For two-body reac- 
tions the rates are computed from Eq.(3). The S- factor 
is a smooth function of energy except for narrow reso- 
nances. When it is assumed to be constant, the integrand 
in Eq.(3) is peaked at a most effective energy which can 
be approximated by a Gaussian function and the inte- 
gral in Eq.(3) can be calculated analytically. In [l2j], 
this approximation is not used and the integral is solved 
numerically for non-resonant as well as broad resonant 
contributions. Also a detailed analysis of uncertainties 
is provided for more realistic lower and upper bounds to 
the adopted rates. The analysis in [ijj was performed 
using polynomial expansions for the cross-sections, and 
the uncertainties on rates were in general only estimated 
by allowing 5-factor limits to encompass all existing data 
whereas in the work of [13] the cross-sections are analyzed 
in the R-matrix framework which provides a more rigor- 
ous energy dependence, based on Coulomb functions. In 
addition to that the advantage of all the available data 
is taken to constrain the S*-factor, not restricting the 
data sets only to the energy range of BBN. Also careful 
evaluation of the uncertainties associated with the cross- 
sections and reaction rates is done based on standard 
statistical techniques. More accurate results are, there- 
fore, expected. The Maxwellian-averaged thermonuclear 
reaction rates of relevance in astrophysical plasmas are 
evaluated [l2[ assuming either that the target nucleus is 
in its ground state, or that the target states are thermally 
populated following a Maxwell-Boltzmann distribution, 
except some cases involving isomeric states. Table-I lists 
nuclear reactions for which Maxwellian-averaged ther- 
monuclear reaction rates are modified [TH, E3] for BBN. 



IV. EFFECT OF THERMONUCLEAR 
REACTION RATES ON PRIMORDIAL 
ABUNDANCES 

A comprehensive study of the effect of nuclear reac- 
tion rates on primordial nucleosynthesis is performed. 
In Fig.l and Fig. 2 the primordial abundances D/H and 
3 'He i 'H of D and 3 He nuclei with respect to the num- 
ber of H nuclei are shown as functions of evolution time. 
The continuous line represents the standard BBN results 
whereas the dashed line represents the same with mod- 
ified Maxwellian-averaged thermonuclear reaction rates 
[l2l ]. In Fig. 3 the 4 He primordial mass fractions are plot- 
ted against the evolution time for standard BBN and with 
the modified reaction rates [l2| • In Fig. 4 the primordial 
abundances 7 Li/ H of 7 Li nuclei with respect to the num- 
ber of H nuclei are shown as functions of evolution time 




Time(s) Time(s) 



FIG. 1: (Color online) Plots of the deuterium abundances ver- FIG. 3: (Color online) Same as Fig.l but for 4 He mass frac- 
sus time. The continuous line represents the standard BBN tions. 
results whereas the dashed line represents the same with mod- 
ified Maxwellian-averaged thermonuclear reaction rates 
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FIG. 5: (Color online) Plots of all the abundances with mod- 
ified Maxwellian-averaged thermonuclear reaction rates [l2| . 
The curve marked 'rest' represents sum of all the abundances 
of nuclei higher than 7 Be. 



for standard BBN and with the modified reaction rates 
[l2j . The plots of calculations incorporating R- matrix 
reaction rates 17] are not shown in the figures as they 
are limited to only small Tg values and do not alter the 
results. In Fig. 5 plots of abundances of all the elements 
with respect to the number of H nuclei are shown as func- 
tions of evolution time. The impact on BBN of the recent 
compilation of thermonuclear reactions rates 



12j does 



not make large overall changes. These modified rates do 
affect the magnitude of these predictions at intermediate 
times for cases such as 4 He or 7 Li. However, as obvi- 
ous from the figures the final values of these predictions 
remain almost same. 



TABLE II: Yields at the test (tj 10 = 3.162) and the WMAP 
(771a = 6.19 ± 0.15) [2l( baryonic densities. 



deduced using the test value (7710 = 3.162) and from the 
WMAP HH results (7710 = 6.19) is provided. In this com- 
parison, the results corresponding to the present work use 
all the thirty- five modified reaction rates [lU, [I?} listed in 
Table-I. We find that using these new reaction rates re- 
sults in insignificant increase in helium mass fraction over 
that obtained previously in standard BBN calculations. 
One observes that even with considerable nuclear physics 
uncertainties, most of these nuclear reactions have mini- 
mal effect on the standard BBN abundance yields of 7 Li. 



V. SUMMARY AND CONCLUSION 

In summary, we find little effect on the standard BBN 
abundance yields by replacing the Maxwellian-averaged 
thermonuclear reaction rates by new ones, even given a 
fair uncertainty in issues that bear on key reaction rates. 
The addition of some new reactions to the BBN code 
had virtually no effect on the BBN abundances [26| . The 
chances of solving either of the lithium problems by con- 
ventional nuclear physics means are unlikely and, if these 
problems stand up to future observations, we may be 
forced into just such non-standard BBN scenarios. At 
present, however, theoretical predictions of the primor- 
dial 6 Li abundance are extremely uncertain due to dif- 
ficulties in both theoretical estimates and experimental 
determinations of the 2 H(a,7) 6 Li radiative capture re- 
action cross section. We also argue that present obser- 
vational capabilities do not yet allow the detection of 
primeval 6 Li in very metal-poor stars of the galactic halo. 
However, if the critical cross section is very high in its 
plausible range and the baryon density is relatively low, 
then improvements in 6 Li detection capabilities may al- 
low the establishment of 6 Li as another product of BBN. 
It is also noted that a primordial 6 Li detection could help 
resolve current concerns about the extragalactic D/H de- 
termination. It could be conjectured, however, that the 
new thermonuclear reaction rates used here may turn out 
to be important for nonstandard BBN scenarios [27l . l28j 
with new particle physics. 





Kawano- 


This work 


Kawano- 


This work 


Observations 


Factor 




Wagoner 




Wagoner 










[Test] 


[Test] 


[WMAP] 


[WMAP] 






4 He 


0.2410 


0.2411 


0.2479 


0.2479 


0.232-0. 258[22] 


xl0° 


D/H 


7.250 


7.277 


2.519 


2.563 


2.82i°;f 9 [23j 


xlO" 5 


3 He/H 


1.546 


1.613 


1.033 


1.058 


0.9-1.3[24] 


xlO" 5 


7 Li/H 


1.268 


1.367 


4.627 


5.019 


1.1±0.1[25J 


xlO- 10 



All these calculations described so far for the standard 
BBN and with the modified reaction rates [l2| were per- 
formed with a test value for the ratio of the baryons to 
photons 77 = 1710 x 1CP 10 = 3.162 x 10~ 10 which repro- 
duces the observed 7 Li abundance for the standard BBN. 
In Table-II, the comparison between BBN abundances 



6 



F. Hoyle and R. J. Tayler, Nature (London) 203, 1108 
(1964). [17 
R. Wagoner, W. A. Fowler, and F. Hoyle, Astrophys. J. 
148, 3 (1967). 

R. Wagoner, Astrophys. J. Supp. 18, 247 (1969). [18 

F. Iocco, G. Mangano, G. Miele, O. Pisanti, and P. D. 
Serpico, Phys. Rep. 472, 1 (2009). [19 
K. M. Nollett and S. Buries, Phys. Rev. D 61, 123505 
(2000). 

R. H. Cyburt, B. D. Fields, and K. A. Olive, Astropart. [20 
Phys. 17, 87 (2002). 

R. H. Cyburt, Phys. Rev. D 70, 023505 (2004). 

R. H. Cyburt, B. D. Fields and K. Olive, J. Cosm. As- [21 

tropart. Phys. 11, 12 (2008). 

P. D. Serpico, S. Esposito, F. Iocco, G. Mangano, G. [22 
Miele, and O. Pisanti, J. Cosmol. Astropart. Phys. 12, 
010 (2004). [23 

G. M. Fuller and C. J. Smith. larXrv: 1009.02771 
G. R. Caughlan, and W. A. Fowler, Atom. Data Nucl. [24 
Data Tables 40, 283 (1988). [25 
C. Angulo et al., Nucl. Phys. A 656, 3 (1999). 
C. J. Copi, D. N. Schramm and M. S. Turner, Science [26 
267, 192 (1995). 

C. J. Copi, D. N. Schramm and M. S. Turner, Phys. Rev. [27 
Lett. 75, 3981 (1995). 
L. Kawano, FERMILAB Report No. PUB-92/04-A, Jan- [28 
uary 1992 (unpublished). 

R. N. Boyd, An Introduction to Nuclear Astrophysics 



(University of Chicago, Chicago, 2008) , 1st ed. 

P. Descouvemont, A. Adahchour, C. Angulo, A. Coc, E. 

Vangioni-Flam, Atom. Data Nucl. Data Tables 88, 203 

(2004). 

M.S. Smith, L.H. Kawano and R.A. Malaney, Astrophys. 
J. Suppl. 85, 219 (1993). 

J. M. Blatt and V. F. Weisskopf, Theoretical Nuclear 
Physics (John Wiley & Sons, New York; Chapman & 
Hall Limited, London.) 

Tapan Mukhopadhyay, Joydev Lahiri and D. N. Basu, 
Phys. Rev. C 82, 044613 (2010); ibid Phys. Rev. C 83, 
039902(E) (2011). 

E. Komatsu et al. [WMAP Collaboration], Astrophys. J. 

Suppl. 192, 18 (2011) [arXiv: 1001. 4538] , 

K. A. Olive and E. Skillman, Astrophys. J. 617, 290 

(2004). 

M. Pettini et al, Mon. Not. R. Astron. Soc. 391, 1499 
(2008). 

T. Bania, R. Rood and D. Balser, Nature 415, 54 (2002). 

A. Hosford, S. G. Ryan, A. E. Garcia-Perez, J. E. Norris 
and K. A. Olive, Astron. Astrophys. 493, 601 (2009). 
R. N. Boyd, C. R. Brune, G. M. Fuller and C. J. Smith, 
Phys. Rev. D 82, 105005 (2010). 

B. D. Fields, S. Dodelson and M. S. Turner, Phys. Rev. 
D 47, 4309 (1993). 

S. Dodelson, G. Gyuk and M. S. Turner, Phys. Rev. D 
49, 5068 (1994). 



